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The structure of neutrons, protons, and other strongly interacting particles is now be- 
ing calculated in full, unquenched lattice QCD with quark masses entering the chiral 
regime. This talk describes selected examples, including the nucleon axial charge, struc- 
ture functions, electromagnetic form factors, the origin of the nucleon spin, the transverse 
structure of the nucleon, and the nucleon to Delta transition form factor. 

1. Introduction 

The challenge in using lattice field theory to solve nonperturbative QCD and un- 
derstand hadron structure from first principles has been to solve full (unquenched) 
QCD in the chiral regime of sufficiently light quark masses that one includes the 
dynamics of the pion cloud and can reliably extrapolate to the physical quark mass. 

During the past year, our collaboration has entered the chiral regime by un- 
dertaking a series of full QCD calculations 1 combining computationally economical 
staggered sea quark configurations generated using the Asqtad improved action by 
the MILC collaboration 2,3 with domain wall valence quarks that have chiral sym- 
metry on a discrete lattice. Use of these configurations has enabled us to treat pion 
masses as light as 359 MeV in volumes with spatial dimension as large as 3.5 fm. 
Initial results are presented below for five pion masses, 359, 498, 605, 696, and 775 
MeV, as well as comparisons from the heavy pion world with SESAM full QCD 
configurations 4 using Wilson quarks at pion masses of 744, 831, and 897 MeV. 

These proceedings will briefly summarize the experimental observables that are 
calculable on the lattice and describe selected recent results. More details may be 
found in recent reviews 5,6 ' 7 , recent publications of our group 8,9 ' 10 ' 11,1,12 and of the 
QCDSF collaboration 13 ' 14 ' 15 ' 16 . 
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2. Nucleon Structure 

Since asymptotic freedom renders QCD corrections to high energy scattering small 
and calculable, high energy lepton scattering provides precise measurements of ma- 
trix elements of the light-cone operator 

where n is a unit vector along the light-cone. Expanding O(x) in local operators 
via the operator product expansion generates the tower of twist-two operators 

whose matrix elements can be calculated in lattice QCD. 

The familiar quark distribution q(x) specifying the probability of finding a quark 
carrying a fraction x of the nucleon's momentum in the light-cone frame and the 
(n — l) th moment of this distribution are specified by the diagonal nucleon matrix 
elements: 

(P\0(x)\P) = q{x), (P|Oj" lW -^ } |P> oc jdxx n ~ 1 q{x). 

Analogous expressions in which the twist-two operators contain an additional 75 
measure moments of the longitudinal spin density, Aq(x). The generalized parton 
distributions H(x, £, t) and E(x, £, t) 17 > 18 > 19 > 20 are measured by off-diagonal matrix 
elements of the light-cone operator 

(P'\0(x)\P) = {(jh))H(x, C, t) + i^((a^n a })E(x, £, t), 

where A^ = P 1 ^ - P»,t = A 2 , £ = -n ■ A/2, and ((r)) = U(P')TU(P) for the 
Dirac spinor U. Off-diagonal matrix elements of the tower of twist-two operators 
(P'|C>q MlM2 '^l-P) yield moments of the generalized parton distributions, which in 
the special case of £ = 0, are 

/ dxx n - 1 H(x,Q,t) = A nfi (t), J dxx n - 1 E{x,Q,t) = B nfi (t), 

where A n ^(t) and B n ,i{t) are referred to as generalized form factors. The physical 
observables considered in this work are special cases of these general expressions. 

Nucleon axial charge 

The nucleon axial charge, gA = (1)a 9 = / dxAq(x) oc f dx(q^^ 5 q), is a fun- 
damental property of the nucleon governing j3 decay. It is an ideal test of lattice 
calculations for several reasons. The isovector combination has no contribu- 

tions from disconnected diagrams, it is accurately measured by neutron (3 decay, and 
the functional dependence on m 2 is known from chiral perturbation theory 21 . Fur- 
thermore, renormalization to account for the difference between regulating quantum 
field theory with a lattice cutoff and in the continuum can be performed accurately 
nonpcrturbativcly using the five dimensional conserved current for domain wall 
fcrmions. Thus, conceptually, it is a "gold plated" test of our ability to calculate 
hadron observables from first principles on the lattice, and in addition, since it is 
known to be particularly sensitive to finite lattice volume effects that reduce the 
contributions of the pion cloud, it is also a good test for an adequate lattice volume. 
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Fig. 1. Nucleoli axial charge, g^- The left plot shows our full QCD calculations in spatial box sizes 
1.5 fm (right three error bars), 2.6 fm (next five points) and 3.5 fm (second error bar superposed 
on the left-most point), compared with experiment (data point at far left). The solid line shows 
a fit to the data using chiral perturbation theory and the dashed lines indicate the error band in 
this fit. For comparison, the right graph adds the results by RBCK in a 1.9 fm box (three lowest 
new points that touch the solid line) and by QCDSF/UKQCD (remaining new points). 



The left panel of Figure 1 shows the results of our calculations in the heavy 
and light pion regimes, together with a fit to a curve defined by chiral perturbation 
theory 21 that enables us to extrapolate the calculations to the physical pion mass. 
Note that the extrapolation goes through the experimental point, and that the 
dashed lines denote the errors in the extrapolation arising from the errors in the 
poorly determined counter-term in the chiral expression for qa- Also note that 
at the lowest pion mass, measurements were made in lattice volumes of 2.6 fm 
and 3.5 fm with statistically indistinguishable results, indicating the absence of 
finite volume corrections at this lattice size. For comparison, the right hand graph 
shows the two other calculations of qa in full QCD. The three points by the RBCK 
collaboration 22 in a 1.9 fm box are consistent with our results within statistics, 
whereas the QCDSF/UKQCD results in 1.5-2.2 fm boxes lie uniformly lower than 
our results and experiment. The axial charge is the most completely analyzed and 
theoretically controlled of our calculations in the light quark regime, and clearly 
demonstrates the quantitative potential of lattice QCD with the emerging level of 
computational resources. 

Quark momentum fraction 

The quark momentum fraction, (x) q — J dxxq(x) cx J dxqj^D^q, is particularly 
interesting physically, because it reflects the fact that a large fraction of the momen- 
tum is carried by gluons rather than quarks. Figure 2 shows the flavor-nonsinglet 
difference between the momentum fraction of up and down quarks, which has no 
contributions from presently uncalculated disconnected diagrams and thus may be 
compared directly with experiment. 

In contrast to the case of the axial charge, chiral perturbation theory yields a 
significant dependence of the momentum fraction on the pion mass. The left panel 
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Fig. 2. The left graph shows a curve with chiral perturbation theory behavior at small pion 
masses adjusted to fit full QCD and quenched results at large pion masses. The right plot shows 
preliminary full QCD calculations extending down to the chiral regime as described in the text. 
The left-most data point, denoted by an open circle, corresponds to a 359 MeV pion mass and 
volume 3.5 fm, and the diamond denoting the same mass in a 2.5 fm volume has been shifted 
slightly to the right for visual clarity. The star at the left in each graph denotes the experimental 
value. 



indicates this dependence by extrapolating calculations in the heavy pion regime 
using the functional form 11 a[l - < ^j0j^yr E l n ( m ™+^i )] + bm l- 

The right panel shows the results with chiral valence quarks and MILC configu- 
rations at the five light masses as well as the heavy quark SESAM results denoted 
by open diamonds. The light quark calculations are renormalized by calculating 
in perturbation theory the ratio of the renormalization factor for the operator of 
interest to the corresponding renormalization factor for the axial current and mul- 
tiplying by the nonperturbative renormalization constant for the axial current 12 . 
Figure 3 shows the ratio of the spin averaged momentum fraction, (x) q , to the spin- 
dependent fraction {x)Aq- In this case, the ratio is nearly mass independent and 
yields excellent agreement with experiment. 

Other observables 

Space limitations preclude showing other results presented in this talk, so we 
shall refer briefly here to two other proceedings where these results have recently 
been published 5,6 . 

The electromagnetic form factors Fx and F 2 , corresponding to A10 and B\q 
defined above, characterize the spatial distribution of charge and current at low 
momentum transfer and the ability of a single quark to absorb a large momentum 
transfer and remain in the ground state. Figure 1 of Ref. 6 shows how the slope 
of Fi increases as the pion mass decreases, reflecting the increase in the spatial 
extent of the pion cloud and how the rms radius extracted from this slope lies 
on a curve parameterized in terms of chiral perturbation theory that is consistent 
with experiment at the physical pion mass. Figure 2 of Ref. 6 shows that the Q 2 
dependence of F2 is consistent with that predicted by the next to leading order 
light-cone wave function 23 F 2 ~ Filog 2 (Q 2 /A 2 )/Q 2 , which agrees with recent Jlab 
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Fig. 3. QCD predictions for the ratio of the spin-averaged momentum fraction, (x) q , to the spin- 
dependent momentum fraction, (x)/\ q , compared with the experimental result (denoted by the 
star at the left). 



data 24 . 

The total quark contribution to the nucleon spin 25 is given by the extrapolation 
to t = of A2Q d (t) and B^it), and the contribution of connected diagrams is 
shown in Figure 5 of Ref. 5 for the case of an 897 MeV pion. Combined with the 
calculation 8 of AS = (1)au + (l)Ad, which yields the contribution of the quark spin 
to the nucleon spin, we learn that considering connected diagram contributions in 
the heavy quark regime, 68% of the nucleon spin arises from quark spin, and 0% 
arises from quark angular momentum. Similar results have been obtained in Refs. 
13,26 _ Evaluation of disconnected diagram contributions and extension to the light 
quark regime will reveal the full origin of the nucleon spin. 

Burkardt has shown 27 that the generalized parton distribution H(x,Q,t) is the 
transverse Fourier transform of the quark distribution, q(x,rj_), so the generalized 
form factor, which can be calculated on the lattice 7,9 ' 10 ' 14 , measures its moments: 

Anfii-Al) = J d 2 r ± J dx x^qix, rj.) e if ^ . 

Physically, we expect that as x — ► 1, the transverse distribution approaches a delta- 
function S(r±), and the slope of the form factor therefore approaches zero. Figure 
6 of Ref. 5 bears out this expectation, showing that the slope of the form factor 
A™~ d (t) decreases dramatically for higher moments weighting large x, and that the 
transverse rms radius decreases strongly with the average value of x. 

The experimental method of choice to reveal the presence of deformation in the 
low-lying baryons is measuring the N - A transition amplitude, where non- vanishing 
electric quadrupole (E2) and Coulomb quadrupole (C2) amplitudes are signatures 
of deformation in the nucleon, Delta, or both. It is convenient to measure the ratio 
of the electric to magnetic form factors, Rem = —Ge2{<1 2 )/Gmi{<1 2 ) and of the 
Coulomb to magnetic form factors, Rsm = H<7|£c , 2(<7 2 )/2mA<?Mi(g 2 )- Figure 4 of 
Ref. 6 shows the results of a new lattice method that for the first time has the 
precision to measure non-vanishing Rem and Rsm ratios 28 ' 29 , and extrapolation 
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of these quenched results in the heavy pion regime to the chiral limit yields results 
qualitatively similar to experiment. 
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